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Abstract 

The thermal and magnetic properties of spin-1 magnetic chain compounds with large single-ion 
and in-plane anisotropies are investigated via the integrable sit (3) model in terms of the quantum 
transfer matrix method and the recently developed high temperature expansion method for exactly 
solved models. It is shown that large single-ion anisotropy may result in a singlet gapped phase in 
the spin-1 chain which is significantly different from the standard Haldane phase. A large in-plane 
anisotropy may destroy the gapped phase. On the other hand, in the vicinity of the critical point 
a weak in-plane anisotropy leads to a different phase transition than the Pokrovsky-Talapov tran- 
sition. The magnetic susceptibility, specific heat and magnetization evaluated from the free energy 
are in excellent agreement with the experimental data for the compounds Ni(C2HgN2)2Ni(CN)4 
and Ni(Ci H 8 N2)2Ni(CN)4-H 2 O. 

PACS numbers: 75.10.Pq,64.40.Cn 
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I. INTRODUCTION 

Haldane's conjecture that spin-S" chains exhibit an energy gap in the lowest magnon 
excitation for 2S even with no significant gap for 25' odd inspired a great deal of experimental 
and theoretical investigation. Rich and novel quantum magnetic effects, including valence- 
bond-solid Haldane phases and dimerized phases , fractional magnetization plateaux 
^| and spin-Peierls transitions 5] have since been found in low-dimensional spin systems. In 
this light, the spin-1 Heisenberg magnets have been extensively studied in Haldane gapped 
materials ly, Ml. The valence-bond-solid ground state and the dimerized state form the 
Haldane phase with an energy gap |2J. The Haldane gap in integer spin chains may close in 
the presence of additional biquadratic terms or in-plane anisotropies. In particular a large 
single-ion anisotropy may result in a singlet ground state 0, |^| which is significantly different 
from the standard Haldane phase. 

The difference between the two gapped phases appears to arise from the ground state 
and excitations. In the Haldane nondegenerate ground state, a single valence bond connects 
each neighbouring pair to form a singlet. An expected excitation comes from breaking down 
the valence bond solid state where a nonmagnetic state Si = at site % is substituted for a 
state Si — 1. In this way a total spin S = 1 excitation causes an energy gap referred to as 
the Haldane gap. Whereas the large-anisotropy-induced gapped phase in the spin-1 chain is 
caused by trivalent orbital splitting. For a large single-ion anisotropy, the singlet can occupy 
all states such that the ground state lies in the nondegenerate gapped phase. The lowest 
excitation arises as the lower component of the doublet is involved in the ground state. This 
excitation results in the energy gap. 

A number of spin-1 magnetic chain compounds have been identified as planar 
Heisenberg magnetic chains with large anisotropy. These include NifC^HsN^NifCN^ 
(abbreviated NENC), Ni(C 11 H 10 N 2 O) 2 Ni(CN) 4 (abbreviated NDPK) [lO, |ll| and 
Ni(Ci H 8 N 2 ) 2 Ni(CN) 4 -H 2 O (abbreviated NBYC) [12]. This kind of system exhibits a non- 
degenerate ground state which can be separated from the lowest excitation. This gapped 
)hase also occurs in some nickel salts with a large zero-field splitting, such as NiSnCkj -6H 2 
13 1, [Ni(C 5 H 5 NO) 6 ](C104)2 3 and Ni(N0 3 ) 2 • 6H 2 H The theoretical study of these 



compounds relies on a molecular field approximation for the Van Vleck equation 16]. To 
first-order Van Vleck approximation, the exchange interaction is neglected. To obtain a 
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good fit to the experimental data an effective crystalline field has to be incorporated. This 
approximation causes uncertainties and discrepancies in fitting the experimental data. Here 
we take a new approach via the theory of integrable models. 

It recently has been demonstrated |l7[ that integrable models can be used to study real 
ladder compounds via the thermodynamic Bethe Ansatz (TBA) and the exact high 
temperature expansion (HTE) method 19, In this paper we present an integrable 

spin-1 chain with additional terms to account for planar single-ion anisotropy and in-plane 
anisotropy. The ground state properties and the thermodynamics of the chains are stud- 
ied via the TBA and HTE. We show that a large planar single-ion anisotropy results in a 
nondegenerate singlet ground state which is significantly different from the Haldane phases 
found in Haldane gapped materials B, Q| • We examine the thermal and magnetic proper- 
ties of the compounds NENC HQ and NBYC Q. Excellent agreement between our 
theoretical results and the experimental data for the magnetic susceptibility, specific heat 
and magnetization confirms that the strong single-ion anisotropy, which is induced by an 
orbital splitting, can dominate the low temperature behaviour of this class of compounds. 
Our exact results for the integrable spin-1 model may provide widespread application in the 
study of thermal and magnetic properties of other real compounds, such as NDPK jiol lri| 



and certain nickel salts |13l Il4l. Il5l Il6l|. 

II. THE INTEGRABLE SPIN-1 MODEL 

In contrast to the standard Heisenberg spin-1 materials, experimental measurements on 
the new spin-1 compound LiVGe 2 06 [21j and the compounds NENC and NBYC Q 
exhibit unexpected behaviour, possibly due to the presence of biquadratic interaction and a 
strong single-ion anisotropy, making it very amenable to our approach. The axial distortion 
of the crystalline field in the compounds NENC and NBYC results from the triplet 3 A2 9 
splitting. Specifically, the triplet orbit splits into a low-lying doublet (d xy , d yz ) and a singlet 
orbital (d xz ) at an energy A C f above the doublet. Inspired by the high temperature magnetic 
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properties of this kind of material, we consider an integrable spin-1 chain with Hamiltonian 

N N 

n = JH + Dj2(s*) 2 + EY,((s*) 2 -(sy) 2 ) 



N 



(1) 



N 

Ho = ^ j^j • Sj+i + (<% • >Sj+i) 2 } 
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23, 



2J, |2£ 



7io is the standard su(3) integrable spin chain, which is well understood 
Here Si denotes the spin-1 operator at site i, N is the number of sites and periodic boundary 
conditions apply. The constants J, D and E denote exchange spin-spin coupling, single-ion 
anisotropy and in-plane anisotropy, respectively. The Bohr magneton is denoted by \xb and 
g is the Lande factor. We consider only antiferromagnetic coupling, i.e. J > and D > 0. 



A. The ground state at zero temperature 

For the sake of simplicity in analyzing the ground state properties at zero temperature, 
we first take E — 0, i.e., no in-plane anisotropy. In this case Hamiltonian which can 
be derived from the su(3) row-to- row quantum transfer matrix with appropriate chemical 
potentials in the fundamental basis, is integrable by the Bethe Ansatz. The energy is given 
by 

Mi 

s = - J £ , ( ih2 i " DN ° ~ ^9H(N + - JV_), (2) 



j=l [Vj j -r ^ 



where the parameters v j satisfy the Bethe equations 
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n V j V i +2 TT V j V l + 1 

(fe) O-i) i 11 



i=l W i ~ U i ~2 1=1 «j 



(fc)_.,(*-l)_i 11 „,(k) _ Jk) 
1=1 

Mk+i(k) (fc+l) 



x 11^ (fc+i) ; i • ( 3 ) 

1=1 W i +2" 

In the above, k = 1,2 and j = 1, and the conventions = vf^ = 0, M 3 = apply. 

N + , Nq, iV_ denote the number of sites with spin S z = 1, 0, —1 in the Bethe eigenstates. In 
he thermodynamic limit, the Bethe ansatz equations (0) admit complex string solutions 
18l | from which the TBA equations can be derived j3 5 27]. Following the standard TBA 



analysis, we find that the ground state in the zero temperature limit is gapped if the single- 
ion anisotropy D>AJ. The singlet ground state is separated from the lowest spin excitation 
by an energy gap A = D — 4 J. This energy gap is decreased by the external magnetic field 
H. At the critical point H c \ = (D — 4J)/ figg, the singlet ground state breaks down. Due to 
the magnon excitation, the magnetization almost linearly increases with the magnetic field. 
Once the magnetic field is increased beyond the second critical point H c2 = (D+4J) / \isg the 
ground state is fully polarized, i.e., in the M = M S plateau region. The magnetization derived 
from the TBA is shown in figure^ We remark that a gapped phase exists only for anisotropy 
values satisfying the 'strong anisotropy' condition D > 4 J. As shown in Ref. for the 

spin ladders, the magnetization in the vicinity of the critical fields H c i and H c2 depends on 
the square root of the field, indicating a Pokrovsky-Talapov transition. In this regime, the 
anisotropy effects overwhelm the contribution from the biquadratic interaction and open a 
gapped phase in the ground state. 

When E ^ 0, the in-plane anisotropy x 2 — y 2 breaks the z 2 symmetry and weakens 
the energy gap. In the presence of the in-plane anisotropy term E, the energies split 
into three levels with respect to the new basis 0o = |O) and <p± = a± \ — 1 )+ 1 1 ), with 
a± = [i^sgH ± a/ (/2 B gH) 2 + E 2 ]/E. In this basis the eigenvalues of the underlying per- 
mutation operator are the same as the eigenvalues using the fundamental basis. The 
model thus remains integrable. We find that if E < D, there is still a gapped phase 
with gap A = D - 4J - ^/(/j B gH) 2 + E 2 for the region H < H&. Here the critical 
field H c i = \J (D — 4J) 2 — E 2 /^Bg- In this gapped phase the ground state is the non- 
degenerate singlet. Subsequently, when H > H c i the state 0_ gets involved in the ground 
state. At the critical point H&, the phase transition is not of the Pokrovsky-Talapov type 
due to the mixture of a doublet state in the 0_ state. The magnetization increases as 
the magnetic field increases. Past the second critical point H c2 = '(D + 4J) 2 — E 2 //j lB g, 
the singlet state is no longer involved in the ground state. The state 0_ fully occupies the 
ground state. As the magnetic field is increased beyond H c2 , the (normalized) magnetization 
M = H/ \/ H 2 + (E/fisg) 2 gradually approaches M s = 1. These novel phase transition may 
be observed from the low temperature magnetization curve, which can be evaluated from the 
TBA equations at T = 0, as per the example in figure d It shows that the gap sensitively 
depends on the single-ion anisotropy and the in-plane anisotropy. These phase transitions 
disappear at high temperatures. In addition, the inflection point at H = \J D 2 — E 2 / fx B g 
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and M = |yl — (E/D) 2 indicates that the probabilities of the components 4>o and 0_ are 
equal. Moreover, if the exchange interaction decreases, the magnetization in the vicinity of 
the critical point H cl increases steeply. For J = 0, i.e. the case of independent spins, the crit- 
ical points H c \ and H c2 merge into one point, at which a discontinuity in the magnetization 
occurs. For D < 4J + E, there is no gapped phase. 



B. Magnetic properties at high temperature 

In order to study thermodynamic properties, we adopt the Quantum- Transfer-Matrix 
(QTM) approach j^q]. Explicitly, following jscj the eigenvalue of the QTM for the model 
(JTJ) (up to a constant) is given by 



Ql(v-%)Q2(V) 

Qi(v-i)Q 2 {v-i) 



+ e^(v)Mv + i) Q ^ (4) 
In the above equation the chemical potential terms are 

A*i = VbqH, \i % = D, fi 3 = -fi B gH, (5) 



where for the moment we take E — 0. We have adopted the notation from [2Cj with 
(f)±(v) = (v±iu N )f, Q a {v) = llttiiv-v^) for a = 1,2, and Qp(u) = 1. Here u N = -J/3/N 
where N is the Trotter number. Following the HTE scheme |2(| , we derive the high temper- 
ature expansion for the free energy of model in powers of J/T. Because the expansion 
parameter J/T is small for weak intrachain coupling J, we may expect the free energy to 
accurately describe the thermodynamic quantities at sufficiently high temperatures, even 
for a small number of terms. To third order, the result is 



i/(r, H) = In Co + C£ 3 - + C* (£j 



Qnl^f • - (6) 
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The coefficients C£ , b = 1, 2, 3, are given by 



^1,0 ~~ 2A + , 

C\ fi = 3A + (l-2A + ) + 3A_, (7) 
C\ fi = fA + (l - %A+ + 8A 2 + ) + 8A_(1 - 3A+), 



with 



Co — Bo t r>, 

A+ = B Dfi /B^ D , 

A_ = eMD/T)/Bl D , 

B x>y = 2exp(x/T) cosh{n B gH/T) +exp(y/T). (8) 

For later use, we also give the HTE free energy with in-plane rhombic anisotropy E. If 
the external magnetic field is parallel to the z-axis, the chemical potentials in Eq. (jlj) become 

Mi = h, fi2 = D, /i 3 = -h, (9) 
where h = \J E 2 + {jj,B9\\H) 2 . In this case the function B xy in Eq. (jHJ) changes to 

B x>y = 2exp(x/T) cosh(/i/T) + exp(y/T). (10) 

On the other hand, if we apply a perpendicular magnetic field to the Hamiltonian ( | T ]) . 
the chemical potential terms in Eq. (j^J) are replaced by 

Mi = \{D-E + ti), 

M2 = E, (11) 
Ms = l(D-E-h'), 



where /i' = y/ (D + E) 2 + Ag\fi 2 B H 2 with a = x or |/ |22|| . Subsequently, we have 



C*0 — -B(D-£)/2,D; 



-4+ — B(d+e)/2,d-e/B(d-E)/2,di (12) 
A_ = exp( J D/T)/ J Bf D _ s)AD , 

with now .B^ = 2exp(x/T) cosh(M'/T) + exp(y/T). 



Eq. (JBJ) for the free energy f(T, H) is our key result. Physical properties such as the 
susceptibility, magnetization and the specific heat follow in the usual way by differentiation. 
We also use f(T, H) to calculate the phase diagram for both T ~ K and finite temperatures 
(see figure EJ). We find that considering up to 3rd order in J/T is sufficient as higher orders 
are negligibly small. This is in stark contrast to other series expansions which need many 
orders to accurately describe physical properties. This is mainly because here the coefficents 
are not just constants, but functions of the external model parameters, e.g., the magnetic 
field and the coupling strength. 



III. SPIN-1 COMPOUNDS 



A. The compound NENC 

It is known that antiferromagnetic spin-1 chains |6j,Uj with weak planar anisotropy can 
exhibit a non-magnetic gapped phase. The large D gapped phase has been observed in the 
compounds NENC, NDPK and NBYC In these compounds the in-plane anisotropy 

x 2 — y 2 breaks the z 2 symmetry and weakens the planar anisotropy. From experimental anal- 
ysis, it was inferred that the in-plane anisotropy E in NENC flOj is negligible in comparison 
with the large D single-ion anisotropy, where the Nickel(II) z 2 orbit along the c-axis forms a 
strong crystalline field. As a result the low temperature physics is dominated by this strong 
crystalline field. The antiferromagnetic exchange interaction further lowers the energy but 
its contribution to the ground state as well as the low-lying excitations is minimal. As a 
consequence, the Hamiltonian (0) can be expected to describe this compound quite well. 
Experimentally, the specific heat was measured up to a temperature around 10 K in the 
absence of magnetic field ^(|. A typical round peak for short range ordering at T « 2.4 K is 
observed, see Figure El An exponential decay is detected for temperatures below approx 2.4 
K. Our calculated HTE specific heat for the Hamiltonian with best visual fit constants 
J = 0.17 K and D = 6.4 K in the case (j3J) (the solid line in FigureEJ) is in excellent agreement 
with the experimental curve in the temperature region T > 0.8 K. In particular, the analytic 
result for the specific heat gives a better fit with experimental data than the result from 
perturbation theory [ljj]. For low temperatures (below 0.8 K), paramagnetic impurities and 
a small rhombic distortion are the main reasons for the discrepancy. The inset of Figure El 
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shows that the inclusion of a small rhombic anisotropy E = 0.7 K gives a better fit for low 
temperatures than with E = 0. However, at high temperature this rhombic anisotropy is 
negligible. 

As far as we know, the susceptibility was measured only for powdered samples of this 
compound. Moreover, the experimental susceptibility of NENC was studied only in the 
temperature range 50 mK - 18 K under a static magnetic field H = 0.1 mT. From the 
data shown in Ref. [10] we cannot accurately estimate the contributions for the Curie- Weiss 
term and the paramagnetic impurity. In Figure HU we present our theoretical curves for the 
susceptibility with parallel and perpendicular field evaluated from the free energy associated 
with different chemical potentials. A susceptibility estimation for powdered samples using 
XPowdcr ~ §X|| + §X± jig] does not fit the experimental data very well at low temperatures 
due to the Curie- Weiss contribution and paramagnetic impurities. A visual fit with the 
experimental susceptibility suggests that the contribution from the Curie- Weiss term is not 
negligible. We find that our theoretical susceptibility XPowdcr for powder with a Curie- Weiss 
contribution c/ (T — 6) gives a satisfactory agreement with the experimental curves, where 
c 0.045 cm 3 K/mol and 8 « —0.9 K. This fit suggests the values J = 0.17 K and D = 6.4 
K, with g± = 2.18 and g» = 2.24. From the TBA analysis we find an energy gap A « 5.72 K 
with a parallel external magnetic field at zero temperature for these coupling constants. The 
typical antiferromagnetic behaviour of the susceptibility with a parallel magnetic field to the 
axis of quantization follows from our results. This is in accordance with the behaviour of the 
specific heat given in Figure El The inset of Figure El shows the magnetization of a powdered 
sample at T = 4.27 K. It is obvious that the singlet is supressed by the temperature. Fitting 
suggests the empirical relation M Powdor « \M\\ + |Mj_ for the powdered magnetization with 
the same constants as before. Here Mm and Mj_ denote the magnetization with the field 
parallel and perpendicular to the axis of quantization. 



B. The compound NBYC 

We now turn to theproperties of the compound NBYC, which has also been experi- 
mentally investigated [l2|. In particular, in-plane anisotropy E and a large anisotropy D 
are present, suggesting that the model Hamiltonian (PQ) may again be a good microscopic 
model for this type of compound. Theoretical studies based on strong-coupling expansion 
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methods |30j suggest that the anisotropy of this compound might lie in the vicinity of the 
boundary between the Haldane and field-induced gapped phases 3]. However, due to the 
validity of the strong-coupling expansion method, the fits for specific heat, susceptibility and 
magnetization become increasingly inconsistent with each other as the rhombic anisotropy 
increases. Figure 0] presents the susceptibility for this compound. The theoretical suscepti- 
bility curve for the powdered sample is evaluated from the free energy © with parallel and 
perpendicular fields (see Eq. (TT2*|) ) via the empirical formula XPowder ~ + §X±- A good 
fit for the susceptibility suggests the values D = 2.62 K, E = 1.49 K, J = 0.35 K, with 
9\\ = g± = 2.05. A small discrepancy at low temperature can be attributed to a Curie- Weiss 
contribution term. From the TBA analysis we conclude that the ground state is gapless. 

The inset of Figure 0] shows the magnetization for powdered samples at 5 K, 10 K and 
20 K. Again our theoretical curves are evaluated using the empirical relation Mp ow d er ~ 
|M|| + |Mj_ for the powdered magnetization. An overall agreement in magnetization for 
different temperatures gives a consistent parameter setting for the susceptibility. The singlet 
state is now supressed by the in-plane rhombic anisotropy and the temperature. 

The specific heat was measured up to a temperature of 6 K in absence of magnetic field 



12j. The theoretical specific heat evaluated from the model Hamiltonian (0) (the solid 



line in Figure EJ), with the same parameters used before, is in good agreement with the 
experimental curve in the temperature region 0.5 K to 6 K. For temperatures below 0.5 K 
the high temperature expansion does not converge and thus cannot provide valid predictions. 



IV. CONCLUSION 



We have investigated the thermal and magnetic properties of spin-1 compounds with large 
single-ion anisotropy, such as NENC and NYBC, via the thermodynamic Bethe Ansatz and 
the high temperature expansion for the integrable model (1). Excellent agreement was found 
with the experimental magnetic properties of these compounds |3l(. The large single-ion 
anisotropy results in a nondegenerate singlet ground state which is different from the valence 
bond solid Haldane phase. The in-plane anisotropy weakens the energy gap. 

Finally, we give the full phase diagram of the compound NENC in Fig. El We see that 
the gapped phase is quickly exhausted as the temperature increases. The magnetic ordered 
Luttinger liquid phase lies between the curves defined by H& and H C 2. The ferromagnetic 
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polarized phase is above the curve. The intersection of the critical curves and the H- 
axis indicates the estimated values H c \ « 3.8 T and H C 2 ~ 4.7 T, which coincide with the 
TBA results at T = K discussed in section 111 Al However, for the case where the in-plane 
anisotropy E ^ 0, the critical behaviour is different from the phase diagram of Fig. El In 
this case the fully-polarized phase appears for H » H c2 because the in-plane anisotropy 
mixes the doublet components \S Z = ±1). We anticipate that the exact results for the 
susceptibility and the magnetization of the powdered samples as well as for the compounds 
with parallel and perpendicular magnetic fields may find widespread use in the study of 
their magnetic properties and for identifying the quantum effects resulting from single-ion 
anisotropy. Our analytic approach via the Hamiltonian (1) may thus describe the thermal 
and magnetic properties of other compounds, such as NDPK NiSnCl 6 ■ 6H 2 jl^ . 

[Ni(C 5 H 5 NO) 6 ](C104)2 and Ni(N0 3 ) 2 ■ 6H 2 Q. 
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FIG. 1: Magnetization versus magnetic field H in units of saturation magnetization for Hamiltonian 
(JTJ) with J = 0.2 K, D = 6 K, E = 3.5 and K, g = 2.0 with parallel magnetic field. The solid 
and dashed lines denote the magnetization derived from the TBA with E = 3.5 K and J? = 0K, 
respectively. The magnetization curve for E = 3.5 K indicates different quantum phase transitions 
in the vicinity of H c \ and than the E = K square root field-dependent critical behaviour in 
the absence of the in-plane anisotropy. 
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Temperature (K) 



FIG. 2: Comparison between theory and experiment [10] for the magnetic specific heat versus 
temperature of the compound NENC. The conversion constant is Chte ~ 8Cexp (J/mol-K). The 
solid line denotes the specific heat evaluated directly from the free energy © with the paramet 
ers J = 0.17 K, D = 6.4 K, g = 2.24 and fiB = 0.672 K/T. The inset shows the low temperature 
specific heat. Clearly the inclusion of in-plane rhombic anisotropy E = 0.7 K (dashed line) gives a 
better fit than without rhombic anisotropy (solid line). 



14 



o 


EXP Powder 


HTE Parallel 




HTE Perpendicular 




HTE Powder 




HTE Powder+CW 


\ 








5 10 15 

Temperature (K) 

n 

FIG. 3: Comparison between theory and experiment [TJj for the susceptibility versus temperature 
of the compound NENC at H = 0.1 mT. The fitting curve (solid line) is obtained via the empirical 
relation XPowder ~ gXll + §X-L together with a Curie- Weiss (CW) contribution. The inset shows 
the comparison between theory and experiment ^(| for the magnetization versus magnetic field of 
NENC at the temperature T = 4.27K. A good fit for both the susceptibility and magnetization 
suggests the coupling constants J = 0.17 K, D = 6.4 K, gj_ = 2.18 and gii = 2.24. The conversion 
constants are xhte ~ 0.8123xEXp(cgs/mol) and Mhte ~ 8.5Mexp (10 3 cgs/mol). 
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FIG. 4: Comparison between theory and experiment [12j for the susceptibility versus temperature 
of the compound NBYC. The conversion constants are the same as for NENC. The solid line is 
the susceptibility for the powdered samples with coupling constants D = 2.62 K, E = 1.49 K and 
J = 0.35 K, with g\\ = gj_ = 2.05. The small discrepancy at low temperature might be attributed 
to a Curie- Weiss contribution. The inset shows the magnetizations for powdered samples at 5 K, 
10 K and 20 K. In each case the theoretical results verify the existence of weak exchange coupling 
and in-plane rhombic anisotropy, with a strong single-ion anisotropy. 
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FIG. 5: Comparison between theory and experiment [lj] for the magnetic specific heat versus 
temperature of the compound NBYC. The conversion constant is Chte ~ IOCexp (J/mol-K). The 
solid line denotes the specific heat at H = 0.1 mT evaluated directly from the free energy (|f)J) with 
the same parameters as in Figure @J 
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FIG. 6: Phase diagram for the compound NENC with parameters J = 0.17 K, D = 6.4 K, g = 2.24 
and /i£> = 0.672 K/T. 
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